A Numerical Model of Crossed Andreev Reflection and Charge Imbalance 
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We present a numerical model of local and nonlocal transport properties in a lateral spin valve structure 
consisting of two magnetic electrodes in contact with a third perpendicular superconducting electrode. By 
considering the transport paths for a single electron incident at the local F/S interface - in terms of probabilities 
of crossed or local Andreev reflection, elastic cotunneling or quasiparticle transport - we show that this leads 
to nonlocal charge imbalance. We compare this model with experimental data from an aluminum-permalloy 
(Al/Py) lateral spin valve geometry device and demonstrate the effectiveness of this simple approach in 
replicating experimental behavior. 

PACS numbers: 74.45.+c,74.40.Gh,74.78.Na,75.75.-c 



I. INTRODUCTION 

Crossed Andreev reflection (CAR) is a charge transfer 
process whereby an electron incident at a normal metal 
or ferromagnet to superconductor junction may enter the 
superconductor at energies less than the gap energy A through 
simultaneous retroreflection of a hole in a spatially separate, 
nonlocal electrode. Previous work by others has demonstrated 
this effect experimentally 1,2 for a nonlocal or lateral spin 
valve geometry consisting of two separate normal state 
electrodes incident with a third superconducting electrode, 
laterally separated on the scale of the BCS coherence length 
£o- This effect has been observed as a negative nonlocal 
voltage V„2 and negative nonlocal differential resistance 
dV„;/dI, where V„; is measured across the second normal 
electrode (defined as the detector electrode) separate to that 
through which current I is applied (the injector) 3 . CAR has 
attracted considerable recent interest due to the potential to 
create solid state quantum entanglement via the splitting of a 
Cooper pair via CAR, demonstrated by recent experimental 
studies^. 

Previous work has shown that other competing processes 
also contribute to the nonlocal effect - elastic cotunneling 
(EC) by which an electron may tunnel nonlocally via an 
intermediate virtual state in the superconductor and nonlocal 
charge imbalance (CI), an effect produced through the 
creation of charge nonequilibrium in the superconducting 
electrode through quasiparticle injection and diffusion 6 . Both 
of these processes have been suggested to fully or partially 
cancel the CAR effect, by producing an equal but opposite 
contribution in nonlocal voltage, or to dominate over CAR 
in the regime near T c or near the critical field. The 
use of ferromagnetic electrodes, magnetized parallel(P) or 
antiparallel(AP), has been suggested and demonstrated as a 
potential means to separate the spin-dependent CAR and EC 
effects 7 

A number of theoretical models have been developed 
to model the lateral spin valve by taking an analytical 
approach - either solutions of the Bogoliubov-de Gennes 
equations, following the method of Blonder, Tinkham and 
Klapwijk(BTK)&i2 for Andreev reflection at a single N/S 
junction, or by solution of the Usadel equations 1 -^. Such 



FIG. 1. SEM image of device. Py magnetic electrodes (1) and (3) 
intersect with an orthogonal Al electrode (2)-(4). Current application 
and voltage detection is performed as in the text 



studies suffer the deficiency of being untested against 
experimental data or being unable to simultaneously fully 
replicate all the nonlocal effects - particularly those due 
to nonlocal charge imbalance and the negative nonlocal 
resistance associated with CAR. In this paper we present a 
simplified means to model these effects by considering the 
possibilities for a single electron incident at the local (injector) 
N/S interface and demonstrate the effectiveness in such an 
approach in terms of replicating real experimental data. 



II. EXPERIMENT 

In order to provide comparison with the model, devices 
were fabricated for experimental measurement by electron 
beam lithography using a Raith 50 system in a standard 
lateral spin valve configuration. These consisted of two Ni:Fe 
80:20 (Py) electrodes with lxl/im and 2x2^m nucleation 
pads and attached nanowire with either a continuous 300nm 
in width or tapering from 600nm to 300nm (over 400nm), 
contacting a perpendicular 300nm width Al electrode, of 
lateral inter-electrode separation 600-900nm (Fig. [T}. 
Patterning was performed using standard PMMA positive 
resist. Deposition of Py was performed by DC magnetron 
sputtering at 34W/2.5mTorr Ar at system base pressure 
3xl0 _9 torr. Overlay of the Al electrode was by EBL overlay 
patterning followed by Al deposition at 50W/2.5mTorr Ar 



at base pressure 3.2xl0 _8 torr, preceded by an in-situ 40s 
Ar + mill at 410V acceleration voltage to clean the Py 
electrode surface of oxide and residual resist. Growth 
thicknesses of the Py and Al electrodes were 15nm and 30nm 
respectively for all devices. Measurement was undertaken in 
an adiabatic demagnetization refrigerator cooled to 600mK 
with current injection and local voltage detection between 
1-2 in Fig. [TJ and nonlocal voltage detection 3-4, using 
both a DC applied current ld c <100^A and nanovoltmeter 
and AC current 7 ac =0.25/iA with DC offset J^ c in order 
to measure differential conductance via a standard lock-in 
amplifier technique. Low AC frequency of 9.99Hz was used 
to minimize inductively generated nonlocal voltage observed. 
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Differential conductance through the local (injector) 
junction was measured using the lock-in method with 
I ac =0.25/zA and DC offset \d c to 50/^A. The resulting 
data can be seen in Figure [3] and exhibited a subgap 
enhancement of conductance consistent with single junction 
Andreev reflection^. The model of Blonder, Tinkham and 
Klapwijk (BTK) and as modified by Strijkers et. al£^ 
was fitted to the experimental data by a standard x-squared 
minimization process to obtain parameters (see figure caption) 
including spin polarization P, interface parameter Z and 
effective temperature T* that were within physical limits. Two 
methods of fitting were used: with a single gap parameter Ai 
and two parameters Ai representing a region of suppressed 
superconductivity and A2, a higher (bulk) value. Although 
the single- A approach was capable of replication of the peak 
at zero and of high ld c behavior, only the two- A approach was 
able to replicate the finite bias minima observed at Id c =10- 
30/iA in the conductance data. This implied the existence 
of a region of suppression of the superconductivity in the 
Al - or a proximity effect induced superconducting region 
in the injector electrode - potentially arising from the effect 
of the adjacent magnetic electrode. The low value of P 
compared to the P=0. 3-0.4 range expected can be attributed 
to the non-point contact nature of the junction and that the 
polarization returned reflects the effective polarization at the 
junction rather than the bulk value for the Py. 



B. Nonlocal Effect 

Examples of nonlocal measurement data taken for a range 
of devices can be seen in Figures |2]|6] Figure |2] shows the 
nonlocal DC voltage V n i measured for a range of devices 
A-D of inter-electrode separation L 150-900nm measured 
as a function of applied DC current through the injector 
ldc- We include a range of properties for these devices 
in Table UJ A negative nonlocal voltage was observed, with 
a general trend of reduction in amplitude with increasing 
inter-electrode separation L characteristic of a CAR effect 
decaying over £ - This effect can also be seen in a plot of 
nonlocal differential resistance dR n i in Figs [4] for two of the 



FIG. 2. Nonlocal DC voltage V n{ vs. I dc at 600mK 
for 4 devices A,B,C,D of inter-electrode separation 
A= 1 50,B=250,C=600,D=900nm 
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FIG. 3. Local normalized differential conductance G/G n for an 
example injector junction, with BTK model fits using a single-A 
(Ai) and by a two-A (Ai and A2) method. Fitting parameters of 
P=0.182, Ai=0.1491meV, A 2 =0.2892meV, Z=0.145 and T*=1.7K 
were used, the position of the Ai and A2 indicated on the figure. 



Device 


Separation L,nm 


Ai(meV) 


A 2 (meV) 


Z 


P 


Device A 


150nm 


0.120 


0.254 


0.051 


0.135 


Device B 


250nm 


0.119 


0.261 


0.120 


0.114 


Device C 


600nm 


0.115 


0.242 


0.051 


0.136 


Device D 


900nm 


0.114 


0.240 


0.024 


0.139 


Device E 


600nm 


0.148 


0.288 


0.145 


0.182 



TABLE I. Summary of device properties. Ai and A2 values are 
those of the injector junction, obtained by fitting to the model of 
BTK. Also included are the spin polarization P, barrier parameter Z 
and inter-electrode separation L. 
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FIG. 4. Nonlocal resistance dR n i measured via the lock-in technique 
at T=600mK as a function of Idc at the injector for Device A of 
inter-electrode separation 250nm and Device D of inter-electrode 
separation 900nm. Measurements are shown for both configurations 
of injector/detector; the reverse (rev) configuration refers to current 
injection through the untapered (300nm width) electrode 



devices A, D demonstrating a characteristic shape of finite 
bias negative minima and zero Idc peak, attributed to EC 
and CAR-dominated transport regimes 1415 . The effect was 
observed to be stable on repetition and highly dependent on 
the injector/detector selection, reversing the choice producing 
a lower dR n i effect likely arising from variation in local 
junction properties of the injector. The absence of a finite 
dRni at high Idc was potentially due to the elimination of 
nonlocal CI effects at the 600mK measurement temperature 
that would otherwise contribute a finite nonlocal voltage. 



C. Temperature Dependence 

Figure [5] shows the effect of changing temperature on dR n i 
vs. Idc for a separate Device E with L=600nm. No effect was 
observed above T c (at T>1.3K in the figure). A reduction in 
both peak height and negative minima depth was observed, 
due to reduction in nonlocal superconducting effects. A 
shift towards ld c =0 of the zero bias minima was observed, 
likely due to reduction in A with increasing T At higher Id c 
sharp inversions in dR n i were observed, the position of which 
reduced in ld c as temperature increased. We attribute these 
features to nonlocal charge imbalance, given the temperature 
and current dependence and similarity in shape to features 
observed by others- notably in the the work by Cadden- 
Zimansky and Chandrasekhar 16 . 

Figure|6]shows the variation in zero bias (ldc=0) peak with 
temperature from 0.45- 1.5K for a separate device. A finite, 
temperature independent effect was observed as T— >Q which 
we attribute to the elimination of nonlocal charge imbalance 
due to the reduction in the charge imbalance length Aq = 
y/Drq, with D the metal diffusion constant and tq the charge 



FIG. 5. Device C nonlocal resistance dR n i as a function of Idc at 
a range of temperatures. No signal was detected at T>1.4K, with 
T c for the device reached at ~1.35K. The asymmetry in ±1 was a 
product of the measurement method 




FIG. 6. Experimental temperature dependence of nonlocal 
differential resistance dR n i for a separate Al/Py Device E of inter- 
electrode separation 900nm, with slightly higher T c ~1.4K 



imbalance time given by: 



T Q = 



4k B T 
7rA(T)' 



(1) 



with n n the inelastic scattering time. 

At low temperatures, the CI effect is minimized, implying 
the existence of other non-CI effects (i.e. - EC or CAR) to 
create the finite peak. As T— >T C , this finite value is reduced 
to zero with a large inversion in dR n i close to T c . No effect 
was observed above T c , as expected if the measurements were 
attributable to nonequilibrium superconducting processes. 
We attribute this inversion effect also to nonlocal charge 
imbalance, arising from the divergent increase in Aq. Such 



an inversion effect has been seen in studies of local and 
nonlocal CI effects 17,17 . Previous works by Beckmann et 
al.— and Kleine et al. 18 have however shown a different form 
for temperature dependence at zero bias, an increase from a 
low or zero nonlocal resistance at low temperature to a finite 
positive peak as T— *-T c , followed by a sharp drop to zero 
above T c . We attribute the difference in this work to be due to 
the presence of an additional finite effect at low temperature, 
observed as a peak at zero bias and likely arising from EC, 
and to a differing negative contribution to nonlocal resistance 
arising from either nonlocal CI or device properties such as 
the effect of a suppressed T c at the interface region. Further 
discussion of this aspect is given in a later section. 



y=y x 



y=0 



xi A 1 x 2 A 



2 x. 




300nm 



600nm 



y=-1000 



III. MODELING 

In order to understand the experimental effects seen here 
and elsewhere in the literature, we consider a numerical model 
of the processes undertaken by a single electron incident at 
the injector. We initially model the incident electron as a 
single charge packet, interacting only with the electric field 
in the wire and undergoing diffusive motion with a mean free 
path Xf in the injector (ferromagnetic) material and Xs in 
the normal (or superconducting) metal. An initial position 
for the electron is defined at ij=-100Xf with the inter-metallic 
interface at y=yo=Q. We define positions x\ and x 2 as the left 
and right-hand edges of the injector electrode and y\ the width 
of the orthogonal electrode, such that for y <0 the electron is 
confined within X\ to x 2 - We also define position j/j as the 
top edge of the normal metal electrode, such that the electron 
is confined to the region y<J/i- 

For the electron dynamics, we set an initial electron 
velocity Vo from a Fermi distribution function at system 
temperature T at a point xo where x\ < xq < x%, j/=-100Af 
and direction is defined by unit velocity vector v such that: 



FIG. 7. Schematic of model configuration for a lateral spin valve 
geometry, with injector (I) and detector (D) electrodes indicated and 
inter-electrode separation L=600nm 



velocity is reset from the distribution function and new 
pseudo-random values set for the direction vector v. 

We calculate the electric field E(x, y) in the structure 
from the derivative of the numerical solution of the Laplace 
equation for the electric potential <fi, including the relative 
resistivities of the metal electrodes and invoking the 
following boundary conditions, where V is the applied bias to 
the injector: 



—V xi <x<x2, y=-lOOXp 
+V X=-Xl im , V >0 (4) 

x < x\ or x > X2, y — or y = y\ 



dx 



j y <0, x — —X\ or x 



X2 



(5) 



Vo 



VO V: 




(2) 



The values of position xq and component of v are taken 
from a pseudo-random number generator. The system is 
allowed to iterate with time step At, with time t=n t At for n t 
timesteps. The electron acceleration a due to the electric field 
in the wire, the velocity v and new position are calculated at 
each time step, giving a time-dependent position s: 
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-100A F 



V X 

VyQ 



1 



t 2 (3) 



dy 



= | x <xi, y=0, y=yi 



(6) 



where -lOOAi? <<0 and \xu m \ the maximum propagation 
extent of the electron in x which is discussed further in the 
following section A. The solution for V=0.1mV, in addition to 
the electron paths taken through the system for V=0.01-lmV, 
can be seen in Fig. [8] and [9] Electron-electron interactions 
are not considered. We use this free electron method as the 
primary model for electron dynamics for all data shown here. 



A. Normal State 



We set the simulated time step At such that the 
distance traveled by the electron on each step |s| = 

Xp, in practice l-2nm so as to 
We define a scattering time tf 



^/lA^P + IA^P « 

reduce computation time. 

where Xf=\s\tf; when t=nTp, with n an integer, the electron 



In the normal state, an electron incident at y=0 is free to 
continue diffusive motion in the perpendicular electrode, with 
step size equal to the new mean free path Xs in the material. 
We define a position \xu m | whereby the electron is considered 
to have fully escaped the injector electrode region. For -xu m 
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FIG. 8. Electric potential <f> numerically calculated in the device at 
V=0. lmV as a function of position x and y, from which a spatially 
varying electric field E(a;, y) = (d(j>/dx, d(f>ldy) could be calculated 



FIG. 9. (Color online) Electron paths in the normal state of the device 
for N=10 electrons and applied injector bias V=0.1mV (red),0.1mV 
(black) and lmV (blue) 



(to the left in Fig. |7j this electron is considered to contribute 
to the local conductance; for +xu, n the nonlocal. An electron 
exiting to the left is modeled to add +1 arbitrary conductance 
unit to the local total 0i and exiting to the right adding +1 to 
the nonlocal total 6 r . For Nt single electrons passed through 
the system, the total conductance is defined as g; = J^n ®i 
and similarly for g r = J2n ®r- 

For the case of inter-electrode separation L>>Ag, in the 
normal state, the charge current flow in +x towards the 
detector electrode should be zero should be zero as there is no 
path to ground, as a result a counter potential acting to drive 
electrons in the +x direction forms and experimentally it is 
this counter potential that is measured with a high impedance 
meter. This acts to drive electrons in -x, producing no net 
charge current flow as observed experimentally (Fig. |5}. 
Experimentally, this potential is detected via measurement 
using a high impedance meter with no current flow into the 
detector electrode. Without taking this 'aspect into account, 
the model would produce an unphysically high value for 
g r , especially as the applied injector bias V— >0. In this 
work we simulate the potential generated at the nonlocal 
electrode required to drive a second electron of equal and 
opposite velocity, obtained from the effect of a calculated 
counter potential (V>0) or from a distribution function (V=0), 
traveling back into the interface region (in the -x direction). 
For summation over N> 1 electrons, this results in no nonlocal 
conductance (g r =0) in the normal state, matching the situation 
observed experimentally (Fig. |5}. 

Figure [9] shows the paths taken by Nt=10 electrons for 
injector bias V=0.01-0.1mV, illustrating the operation of the 
model in the normal state. For low bias V=0.01mV, the 
electron may take a diffusive path in ±x, indicated by the 
tracks towards +x (right). For high bias, the majority of 
electrons move in -x, the paths confined closer to the inner left 



Process 


8i , Local 


8 r , Nonlocal 


AR 


+2 





CAR 


+1 


+1 


EC 


+1 


-1 


Normal reflection 


-1 





Quasiparticle(CI) 


+1 


±A q c- x °/ A Q 



TABLE II. Single-electron contributions from each process to local 
and nonlocal conductance with zero spin polarization P m . At the 
nonlocal electrode CAR is defined to make a positive contribution 
(exhibited as a negative nonlocal resistance) and EC a negative one 



edge of the junction in the -x direction as they approach y=Q. 
For V=0, electron diffusion can occur with equal probability 
in either direction, zero net charge current in the +x electrode 
imposed as a result of modeled injection of a counter electron 
moving in -x (not shown) representing an opposing counter 
diffusion and ensuring g r =0. However, for V=0 no net flow 
of electrons should occur to the injection region (i.e. Nt=0 
at V=0), which would limit the validity of the model to finite 
bias in the normal state. This was accounted for in the model 
by fixing a finite maximum number of iterative time steps to 
<2xl0 5 , such that the electron never reaches Xu m at V=0 in 
the normal state, giving g r =gi ~0. 



B. Superconducting State 

We next consider the superconducting state by introducing 
two spatially variant values of the superconducting gap 
A\(x,y) and Aa(a;, y) representing a suppressed and bulk 
energy gap, such that A 2 (x,y) > Ai(x,y). This follows 
the approach of Strijkers et al. 13 considering a proximitized 
layer of lower- A at the normal-superconductor interface. For 



the model, Ai=A2=0 for y<0 in the normal state electrodes. 
We restrict the spatial extent of the suppressed Ax to around 
the region of overlap between the normal electrodes and 
superconducting electrode for 0<y<yi and \x^\ with A2 
covering the remainder of the superconducting region, such 
that an electron may progress as a quasiparticle excitation 
in the Ai region for Ai<E<A2 but not into the bulk A2 
region. This is illustrated in Fig. [7] Each electron is 
defined as possessing an energy E, taken at random from 
a Fermi distribution f(E-V) = l/(e^ E - v ^ k » T +l) at each 
scattering step, where f p is a pseudo-random value (0..1), V 
the simulated potential across the injector and T the simulated 
temperature. 

For an electron diffusively reaching the region of nonzero 
A and E<Ai, only transmission by Andreev reflection is 
possible. The energy and potential dependent probability 
for AR Par is calculated from expressions given by from 
the model of BTK 8 and the modification by Strijkers. We 
consider that a fraction of this Andreev reflection may occur as 
CAR such that Pcar=cPar. and that a fraction of the current 
which is not transmitted by AR could be by EC such that 
PEC=Ct{i-PAR+PcAR) defining a probability conservation 
relation P ar+Pc ar.+P n r+P bc=1 including the probability 
of normal reflection Pnr- We assume a similar energy 
dependence for Par and Par. and that the factors c and c t 
defining the CAR and EC fractions are invariant with electron 
energy or applied simulated potential V. For each modeled 
electron, the process undertaken is selected using a pseudo- 
random value p3=(0..1), the range to 1 subdivided according 
to the ratios of probabilities for each process. An electron 
undertaking any of the processes converting to supercurrent 
is said to have exited the system, in a manner similar to that 
of diffusive exit past \xu m \ and it's further diffusive motion 
is not considered - the contributions #,./; for each process are 
given in Table ITU 

For Ai<E<A2 the electron may still Andreev reflect, 
cotunnel or may now enter the superconducting electrode as 
a quasiparticle. We consider the quasiparticle state to be 
electron-like and modeled by random walk diffusive motion 
- but limited to the Ai region within \x%\, with quasiparticle 
transport into the bulk A2 region outside prohibited. For 
the calculation of the electron dynamics, we assume the 
electrostatic force on the quasielectron F(x,y)=Q within 
the superconducting regions, permitting motion with equal 
likelihood in all directions. Such a method for quasiparticle 
diffusion has been previously used in the literature 19 . 

For E>A2 the restriction to \x$\ is removed, permitting 
quasiparticle diffusion throughout the superconducting 
electrode (whilst prohibiting any Andreev reflection). A 
quasiparticle reaching \xu m \ contributes to both the local 
and nonlocal conductance according to Table [II] with an 
exponential decay depending on the separation between 
the point of exit and the nonlocal electrode x s , and decay 
scale Aq the charge imbalance length. This represents 
a phenomenological replication of the nonlocal potential 
produced by charge imbalance, of either negative or positive 
contribution to the nonlocal effect. The same dependence 
on inter-electrode separation is introduced for EC and 



CAR by reducing their fraction of the AR and non-AR 
probability at the injector according to an exponential form 
c=Coe _3:s /'» and c^Ctoe^ 3 ^ on a decay scale equal to 
the BCS coherence length £o> c o and c t are defined at zero 
separation (i.e. - at the injector). Since E is taken from a 
Fermi distribution, even for low V all processes are possible 
at finite T The spatial variation in A permits the simulation 
of regions of energy gap lower than the bulk value for the 
superconducting electrode. 



C. Spin Dependence 

Finally, we consider the case for a ferromagnetic injector 
and detector electrode by considering a total pool of available 
electrons of each spin type: Ni-j-, N14, at the injector and N21-, 
N24,, such that the total number of simulated injected electrons 
Nt=Ni^+N]4 and defining the polarization for the spins in the 
model P m at the injector electrode: 



p — 



(JVit - N u ) 
(N lt + N u ) 



(7) 



P m is assumed to be uniform for both injector and detector and 
represents a means of replicating the bulk spin polarization 
P in the ferromagnetic electrode. A simulated electron is 
selected at random from one of the spin type pools (Ni^/j,)). 
For P m =0 this has no effect on the conductance totals or 
probabilities for any effect. For finite P m , such that Ni-|- ^Ni^ 
and N2-1- 7^24., one pool will be exhausted before the other, 
prohibiting the spin dependent processes (such as AR, CAR) 
from occurring. For example, in the half metallic case where 
Ni-|-=N2f=0, CAR would be prohibited due to the absence 
of the pairing spin-f required to form a Cooper pair. This 
method also permits simulation of magnetization dependence 
by introducing an imbalance in the spin dependent pools such 
that Ni-j- >N]4 and N2f <N2^ in the antiparallel configuration 
and the reverse for parallel magnetization. The model neglects 
spin accumulation and the non-local spin valve effect by 
choice, although such an effect could in principle be included. 



D. Model Examples 



Figures [T0lfT2l show the effects on the modeled nonlocal 
differential resistance g^ 1 for an arbitrary selection of realistic 
parameters (£o=600nm, AQ=l/j,m, Ai=0.2meV, A2=0.4meV, 
T=600mK, T C =1.4K) as a function of inter electrode 
separation and of temperature for Nt=10000 electrons, both 
following qualitative trends observed by others 18 ' 15 . Figure 
[TT1 can be compared directly to Figure |5] from experiment, 
the same qualitative trend of reduction in peak and negative 
minima observed from the model. 

Figure Q~2] shows the temperature dependence at V=0 
(ldc=0) for 3 configurations - no CI (A q =0) and finite 
CAR/EC c=c t =0.2 and with CI switched on: A q =l 
with absence/presence of CAR/EC. We consider the case 
Ai=A 2 =0.2meV with a single T c of 1.4K. We model CI 





FIG. 10. Example of model output I: nonlocal differential resistance 
g^ 1 as a function of V for inter-electrode separation L=300-900nm. 
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FIG. 11. Example of model output II: differential resistance g^T 
vs. V at a range of temperatures T=400mK-4.2K, with A q =0 and 
parameters as in the main text 



as a negative contribution to the nonlocal resistance g^ 1 to 
assist replication of the effects observed experimentally. In 
the case of Cl-only, no finite effect is seen at low-T. This did 
not match the result for the real Al device in Fig. |6] which 
qualitatively follows the model form with finite-CAR/EC 
contribution as T— K) and CI divergence as T— s-T c . Divergent 
behavior is replicated as T->T C in the cases where A q > 
and a finite CI contribution is considered. Assumption of 
a negative CI contribution (A g =l) combined with a finite 
positive low temperature CAR/EC contribution produces a 
negative divergence towards T c but of a form not reflected 
in experimental data. We also note that assumption of a 
positive contribution to the nonlocal resistance (negative A q ) 
with increasing T, the behavior for A 9 =-l, c=c t =0 (CI only, 
dashed line in Figure [T2i i replicates the observed behavior 



FIG. 12. Example of model output III: T-dependence for 
absence/presence of charge imbalance A q =0 or ±1 and 
absence/presence of EC/CAR (c,c t values) at V=0. The result 
for negative CI contribution (-A q , solid) is shown with the mirrored 
result for positive contribution (+A 9 , dashed) 



seen by othersi^^ii as an increase to a finite peak from low-T 
zero g~ l with zero effect above T c . 



Figure Qj] shows the modeled case whereby the 
T C =T^=1.4K of the Ai(T=0)=0.213meV region is set 
0.04K lower than the bulk Tf =1 .44K of the A 2 (0)=0.22meV 
region, such that for T^<T<T,f a fully normal state region 
is formed at the overlap between the ferromagnetic and 
superconducting electrodes. This would be anticipated 
experimentally, especially in the case of Ohmic contacts 
such as utilized in devices studied in this work, due to the 
suppression of superconductivity at the F/S interface. Within 
this temperature range, model nonlocal processes such as 
CAR and EC are limited due to the low gap energy in the 
high-T regime. At low bias within the interface region, 
diffusive electron motion permits incidence with the A 2 
interface at either -x% or +X3. An electron arriving at the 
+23 interface may enter the superconductor primarily as a 
quasiparticle excitation close to T c (due to the low value 
of A2) and with energy E>A 2 - we define this electron 
as producing a negative nonlocal contribution to g^ 1 of 
magnitude \A q \e~ Xs / AQ . In contrast, for the fully normal 
state Ai=A2=0 CI cannot occur to produce a nonequilibrium 
population of electrons and a detectable potential difference, 
as discussed in section A. A combination of the reduction 
in x s , a higher number of electrons reaching +X3 and the 
proximity to bulk Tf producing a longer Aq and larger 
CI effect results in a negative effect that can be far greater 
than that produced by nonlocal CI. This provides a potential 
explanation for the divergent features observed in experiment 
as plotted in Fig. [6] with which excellent qualitative 
agreement is obtained. 
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FIG. 13. Example of model output IV: (solid) the case where the 
region -xi<x<X2, y>0 enters the normal state at T=T,f-0.04K 
with T,f =1.4K the bulk critical temperature of the A2 region which 
remains superconducting. In this case, transport may occur by direct 
quasiparticle transmission (or minimal AR) at the A2 interface in 
±ir. The result is a sharp, negative minima qualitatively matching 
the effect observed experimentally (dashed, scaled from Fig. [6]l 



E. Match to Experiment 

In order to demonstrate the capability of the model, the 
transport behavior of the real device was compared to the 
model result with input parameters realistic for the fabricated 
devices. The result of this is shown in Figure [l4l giving the 
model result overlaid onto real nonlocal differential resistance 
data from experiment. Model parameters were selected to 
best represent the device or from measurement of device 
properties. Parameters xi=-150nm, X2=+I50nm (giving 
width=300nm) were fixed based on the physical dimensions 
of the device. The point of diffusive exit from the system 
|#Zim|=450nm, taken as the half-way point to the detector 
electrode with inter-electrode separation L=600nm, based on 
the assumption arrived at through simulation that the majority 
of electrons crossing past x=300nm would follow a path 
reaching the second electrode at x=750nm. \xu m \ was also 
minimized to reduce computation time, but to be sufficiently 
high to inhibit premature exit of the electron from the system. 

Values of A s =15nm, Ap=5nm and AQ=l/!m were used, 
based on physically reasonable parameters from previous 
work for nonlocal spin value devices and on nonlocal 
charge imbalance. £o=500nm was taken, assuming a 
relatively clean superconducting Al electrode, but with a 
degree of suppression of the coherence length from the 
BCS value and sufficiently long to produce non-negligible 
superconducting nonlocal effects over the inter-electrode 
separation used. Although a value was set for Aq, based 
on the temperature independence of the nonlocal effect at the 
600mK experimental measurement temperature (Fig. [6} the 
assumption was made of complete suppression of nonlocal 
charge imbalance. We explicitly introduce the absence of 
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FIG. 14. Nonlocal differential resistance as a function of applied 
injector DC current Idc for modeled and experimental data. 
Device inter-electrode separation was L=600nm and experimental 
measurement temperature 600mK. 



CI by setting A 9 =0, in order to demonstrate that the effects 
observed can originate solely from non-CI processes. The 
assumption that CI can be completely eliminated in this 
manner was based on the hypothesis from experimental data 
that the origin of the nonlocal effect at T<<T C was from 
processes operating over a decay scale £0 (i.e. - CAR/EC) 
rather than CI. This assumption would not be valid for T— >T C , 
where additional CI effects would arise from the divergence 
of Aq and cannot replicate the temperature dependent effect 
observed experimentally. 

Modeled T c was set to 1.4K, close to that observed 
in experimental data. Using this value of T c , we 
calculated physically reasonable values for A 1 =0.112meV, 
A2(0)=0.301meV derived from standard BCS theory and fits 
to local (single junction) conductance data. We assume the 
simplest case with uniform values of Ai and A 2 for injector 
and detector junctions to the superconductor. The requirement 
for the two values of A was based on the single junction 
experimental data, whereby a region of lower energy gap 
was required to correctly fit the conductance behavior - this 
region likely being that immediately adjacent to the magnetic 
electrodes. On this basis we limit the Ai region extent \x^\ 
to within the overlap region of the injector/detector electrodes 
-150<a;<150nm . We neglect any proximity effect into the 
injector, confining the superconducting region of nonzero A 
strictly to within the perpendicular lateral electrode for 300< 
y <0nm. 

The values of c and c t represented free parameters which 
were selected to give the best fit to the experimental data, each 
parameter effectively controlling negative minima depth and 
zero bias peak height of nonlocal resistance g^ 1 respectively. 
We make the assumption that both c and c t are invariant with 
other parameters such as injector bias. Values of c=0.28 and 
c t =0. 15 were found to give a best fit, implying the presence of 
both effects but an imbalance in a 2: 1 ratio between them such 



as not to cancel each other, with each producing a different 
effect at varying simulated bias. 

In order to translate the arbitrary conductance units used in 
the model to real nonlocal resistance (in Ohms) the conversion 
g r =kSR n i was used with k a scaling factor of approximately 
1000. A total of N T = 10000 electrons were simulated for 
each step in simulated injector V. A good fit was found to the 
experimental conductance data. As can be seen in the figure, 
the only discrepancies arose in matching the height of the 
peak, potentially due to a slightly high values of c or Cj, and 
a discontinuity around 0.23-0.3mV due to the 2-A approach 
taken, rather than a continuous range of A arising from the 
proximity effect. The fit was obtained with nonlocal charge 
imbalance explicitly removed (A 9 =0), in order to demonstrate 
the absence of such effects at low temperature. 



IV. CONCLUSION 

In conclusion, we have fabricated Al/Py devices suitable 
for measurement of nonlocal superconducting effects. We 
observe behavior characteristic of EC or CAR in the low 
temperature regime and nonlocal charge imbalance, arising 



as divergence in nonlocal differential resistance, close to 
T c . In order to further investigate there results we have 
developed a simple numerical model of nonlocal processes 
in a lateral spin valve geometry device that is capable of 
qualitatively replicating the nonlocal effect trends with T 
and bias in experimental work here and in the literature, 
associated with the superconducting processes. The model 
is capable of quantitatively matching the nonlocal resistance 
from a real Al/Py device, using realistic parameters for the 
model variables. The single electron approach offers potential 
extension for shot noise calculation, of interest for recent 
cross correlation and entanglement studies. Although possibly 
inferior to a true analytical model, the numerical approach 
encompasses the key aspects of the physics of such a device 
and gives insight and the ability to replicate the typical local 
properties of a device and provides a simple way to derive 
important parameters such as £ an d Aq. 
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